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AN ALGORITHM FOR COMPUTING THE MULTIGRADED HILBERT 

DEPTH OF A MODULE 

BOGDAN ICHIM AND ANDREI ZAROJANU 

Abstract. A method for computing the muhigraded Hilbert depth of a module was 



^ ' presented in I.14J . In this paper we improve the method and we introduce an effective 

^~^ algorithm for performing the computations. In a particular case, the algorithm may also 



be easily adapted for computing the Stanley depth of the module. Thus, we completely 
^O ■ solve an open problem proposed by Herzog in [.lOJ . 



u 

<[^ ■ 1. Introduction 

^ '. In this paper we introduce an algorithm for computing the Hilbert depth of a finitely 

c^ I generated multigraded module M over the standard multigraded polynomial ring R = 

K[Xi,.. . ,X„]. The algorithm is based on the method presented in [|T4l and some extra 
improvements. It may also be adapted for computing the Stanley depth of M if dim^ Ma < 
1 for all a E Z". As a consequence, we give a complete answer to the following open 
problem proposed by Herzog in [|TOl : 

*vq ■ Problem 1. fTOl Problem 1 .66] Find an algorithm to compute the Stanley depth for finitely 

t^ • generated multigraded R-modules M with diniK Ma < 1 for all a G Z". 

(^ I In recent years, Stanley decompositions of multigraded modules over R have been dis- 



m I cussed intensively. Such decompositions, introduced by Stanley in [I20II . break the module 

M into a direct sum of Stanley spaces, each being of type mS where m is a homogeneous 
element of M, S = K[Xi^ , . . . ,X,^] is a polynomial subalgebra of i? and Sf]Annm = 0. One 

^ I says that M has Stanley depth s, sdepthM = s, if one can find a Stanley decomposition in 

which d > s for each polynomial subalgebra involved, but none with s replaced by 5 + 1 . 

The computation of the Stanley depth is not an easy task, due mainly to its combina- 
torial nature. A first step was done by Herzog, Vladoiu and Zheng in [fT3l , where they 
introduced a method for computing the Stanley depth of a factor of a monomial ideal 
which was later developed into an effective algorithm by Rinaldo in (TW\. Some remark- 
able results in the study of the Stanley depth in the multigraded case were also presented 
by Apel (see ffl, [H), Herzog et al. (see [HI, ^) and Popescu et al. (see [0, liTTll ). 

Hilbert series are the most important numerical invariants of finitely generated graded 
and multigraded modules over R and they form the bridge from commutative algebra to its 
combinatorial applications (we refer here to classical results of Hilbert, Serre, Ehrhart and 
Stanley, see [4J). A new type of decompositions for multigraded modules M depending 
only on the Hilbert series of M was introduced by Bruns, Uliczka and Krattenthaler in 
im and called Hilbert decompositions. They are a weaker type of decompositions not 
requiring the summands to be submodules of M, but only vector subspaces isomorphic 
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2 BOGDAN ICHIM AND ANDREI ZAROJANU 

to polynomial subrings. The notion of Hilbert depth hdepthM is defined accordingly. 
Several results concerning both the graded and multigraded cases were presented in [|9l, 
|fT6]| and ETTl . All of them are based on both combinatorial and algebraic techniques. 

The contain of the paper is organized as follows. In Section |2] we recall some results 
conceming Hilbert depth that will be used in this paper. 

Section [3] is devoted to improve the method presented in [|T4l by restricting as much as 
possible the search for a suited Hilbert decomposition. Theorem [TOl shows the existence 
of upper-discrete Hilbert partitions of degree s for hdepthM > s, so for computing the 
Hilbert depth it may be better to consider only this type of partitions. It is a generalization 
of both ifTSl Lemma 3.4] and [fT9l Lemma 3.3] (notice that, in the particular case of a factor 
of a monomial ideal, the Hilbert partitions coincide with the poset partitions considered 
by Rinaldo and Shen). 

In Section |4] we introduce a recursive algorithm for computing the multigraded Hilbert 
depth of a module (see Algorithm [T]). The algorithm is relative easy to implement because 
of its recursive form and may also be used directly for computing the Stanley depth in the 
case of a factor of a monomial ideal. A non-recursive algorithm for computing Stanley 
depth in the case of a factor of a monomial ideal was introduced in [18, Algorithm 1]. 

Hilbert decompositions are intimately related to Stanley decompositions: All Stanley 
decompositions are Hilbert decompositions; moreover, the latter are prerequisites to the 
existence of Stanley decompositions. In Section [5] we investigate how strong is this con- 
nection. We assume that dimj^M^ < 1 for all a E Z" and we show that Algorithm [T] may 
be easily modified for computing Stanley depth in this case (see Algorithm |2l). This solves 
completely Problem [1] 

We end this section with a vague remark. In the particular case of a normal affine 
monoid, suited Hilbert decompositions have been used with success in order to design 
arguable the fastest available algorithms for computing Hilbert series (see [|51, [|6l and 
lITl). It is an interesting open problem if it is possible to use suited Hilbert decompositions 
in order to design efficient algorithms for computing Hilbert series in other cases. This 
paper is likely to be of particular interest to developers of computer algebra systems. 

2. Prerequisites 

Let R = K[Xi, ...,Xn], with K a field, and let M be a finitely generated Z" -graded R- 
module. In [|T4l the authors presented a method for computing the multigraded Hilbert 
depth of M by considering Hilbert partitions of its Hilbert series. In this section we recall 
the method of [[HI. 

A natural partial order on Z" is defined as follows: Given a,b E Z", we say that a ^b 
if and only if a,- < bt for i = 1 , . . . , n. Note that Z" with this partial order is a distributive 
lattice with meet aAb and join a\/ b being the componentwise minimum and maximum, 
respectively. We set the interval between a and b to be 

[a,b] = {cEZ" \a^c^b}. 

We first recall a definition and a result of Ezra Miller (see [,15,1 ') which will be useful 
in the sequel. Let g eN". The module M is said to be N"-graded if M^ = for a ^N" 
and M is said to be positively g-determined if it is N"-graded and the multiplication map 
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■Xi : Ma — y Ma+ei is an isomorphism whenever a,; > gi. A characterization of positively 
^-determined modules is given by the following. 

Proposition 2. ifTSl Proposition 2.5] The module M is positively g-determined if and only 
if the multigraded Betti numbers ofM satisfy Poa = j8i,a = unless ^ a ^ g. 

Let 

R{-a)Pu, ^ R(^-afo.a ^M^O, 

aeZ" aeZ" 

be a minimal multigraded free presentation of M and assume for simplicity, and without 
loss of generality, that all jSq^a = (and a fortiori all /3i « = 0) if a ^ N". 

Let ^ G N" be such that the multigraded Betti numbers of M satisfy the equalities /3o,a = 
Pi,a = unless ^ a ^ g. Then, according to Proposition |2l the module M is positively 
^-determined. Let 

HMiX)= £ H{M,a)X" 

aeN" 

be the Hilbert series of M and consider the polynomial 

HM{X)^g:= £ H{M,a)X''. 

For a,b eTj" such that a ^b,we set 

Q[aMiX):= £ X' 

and call it the polynomial induced by the interval [a,b]. 

Definition 3. We define a Hilbert partition of the polynomial HM{X)^g to be an expres- 
sion 

^■.HM{X)^8 = Y.QWMiX) 

as a finite sum of polynomials induced by the intervals [a', Z?'] . 

Further we need the following notations. For a ^ ^ we set Z^ = {Xj \ aj = gj}- More- 
over we denote by K[Za] the subalgebra generated by the subset of the indeterminate s Z^. 
We also define the map 

p:{0^a^g}^N, p{a):=\Zal 

and for ^ a ^ Z? ^ ^ we set 

^[a, b] = {ce [a, b] \ Cj = aj for all j eN with Xj e Z^}. 

The main result of [|T4l is: 

Theorem 4. [|T4l Theorem 3.3] The following statements hold: 

(1) Let^: HM{X)^g = E-^j Q[a\b']{X) be a Hilbert partition ofHM{X)^g. Then 



D(q3):M-0( K[Z,,]i-c)) W 



J=l ce«ffa'\fc'l 



4 BOGDAN ICHIM AND ANDREI ZAROJANU 

is a Hubert decomposition ofM. Moreover, 

hdepth D (q3) = min{p (Z?' ) : / = 1 , . . . , r } . 

(2) Let T) be a Hilbert decomposition ofM. Then there exists a Hilbert partition ''^ of 
HM{X)^g such that 

hdepth 2) (<P) > hdepth S). 

In particular, hdepthM can be computed as the maximum of the numbers hdepth 2) (*p), 
where *P runs over the finitely many Hilbert partitions ofHM{X)^g. 

We see that, in order to effectively compute the Hilbert depth of M, we may use the 
following corollary. 

Corollary 5. |[T4l Corollary 3.4] There exists a Hilbert partition 

^■.HMiX)^g=Y^Q[a\b']{X) 
1=1 

of HM{X)^g such that hdepthM = mm{p{b') : /= l,...,r}. 

3. Restricting the search for a good partition 

As seen in the previous section the Hilbert depth ofM can be computed by considering 
all Hilbert partitions of HM{X)^g. In practice, the number of possible partitions can 
easily become huge. For many practical purposes (for example for implementation of the 
method in a computer program) one needs to restrict (as much as possible) the search for 
a partition which will finally provide the right Hilbert depth. In this section we show that 
an improvement is indeed possible. Our results are extending some of the ideas presented 
by Giancarlo Rinaldo in [flSH and Shen in [19] for computations of Stanley depth in the 
case of a factor of a monomial ideal to the general case of a finitely generated module. 

Since many results in this section depend on a number g eN" such that M is positively 
^-determined, we shall assume that g is fixed and known from previous computations (for 
example by using Proposition^. 

Definition 6. Let 5 be a subset of N" and < s <n. We define two subsets of 5, 

B<s'={a^B:p{a)<s} and B>s := {a E B : p{a) > s}, 

where p is the function defined in Section [2l 

Our purpose is to test if M has a partition %i whose hdepth is equal to s. To reach 
this goal set B = {a : X"^ is a monomial of the polynomial HM{X)^g} and consider B as 
a disjoint union of the two sets defined above 

B = 5<i U5>5. 

It is easy to observe that if ^ is a Hilbert partition of HM{X)^g, then we may write 
<P=A+A', suchthat 

A = l^Q[a\b']{X), A'=l^Q[aJ,b^]{X) 

ia ye/' 
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where a' G 5<^ and a^ E B>s for all / G / and j G /'. Then *p can further be refined to a 
new partition ^' = A +A" with 

A"='£Q[aJ,aJ]{X) 
jei" 

where a^ G 5>v for all j G /". 

Therefore, if a partition ^ with hdepth = s exists, then the part A of *p is composed 
of intervals 2[a, b] {X) where a G B^^ and b G B>s- At a first look, in order to find A, we 
have to consider for each element a G 5<s all possible candidates b G 5>.v with a ^ Z?. In 
the following we show that the list of candidates can be reduced considerably. 

Proposition 7. Let P = Q[a, b] {X) be a polynomial such that b ^ g and p{a) < s < p (b). 
Then for each 

b^ G Min{;c : a ^x ^b, p{x) > s} 
there exists a disjoint decomposition ofP 

P = Po + tp,, (*) 

such that Pq is the polynomial induced by the interval [a, b% Pi is the polynomial induced 
by the interval [a\ b% ¥ = b and p{b') > sfor all i = 1 , . . . , r. 

Proof. We see that 

P=(Xf'+...+xf')---(X„"" + ...+X^) 

= X"(l + ... +Xfi-"i) ■ ■■(! + ... +X,^-"«), 

so we may assume for simplicity and without loss of generality that a = (0, ...,0) G N". 
Then we have 





P=(l+Xi + ...+xfi)-- 


■(l+X„ + ...+X* 




^Po+j;^p^, 




where we set 








Po = (l + ...+xf^)-- 


■(i+...+x5) 


and 






Pi = {\ + ...+xf)- 


■■(i+...+x;-)(xf'+^+, 


...+xf')(i + -+. 



*'+n n I I vK 



for all / = 1, . . . ,r (in case Z?^ = b, the term Pi is simply 0). We deduce the following 
description 



M. = ^ ^"' if^'<'' 



•^ 1 bj, otherwise. 
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Since b^ ^b' ^b ^ g,we get that p{b') > p{b^) > 5 as needed. 

We claim that (*) is a partition of [0,b]. To prove this it is enough to show Mon(/',) fl 

r 

Mon{Pj) 7^ if and only if i = j and that the equality P = Po+ £ P, holds. 

r 

For the equality we will show that Mon(/') = Mon(/'o) U |J Mon(P/). We only need to 

(=1 

r 

show that Mon(/') C Mon(/'o) U |J Mon(i',) because the other equality is obvious. 

Let u G Mon(i'), u=X^'. If ci > b^ then u G Mon(/'i), otherwise for sure u ^ Mon(Pi). 
If ci < Z?i , we check if C2 > b*^. If so w G Mon(/'2), otherwise m ^ Mon(i'i) UMon(/'2)- So 
after checking for all variables, we find that u G Mon{Pi) if c/ < b'j, for all 7 = 1, . . . , / — 

1 and Ci > b^ ox u E Mon(Po) if ci < b^ for all / = 1, . . . ,n. It is also clear from this 
description that Mon(/',) fl Mon{Pj) 7^ if and only if / = j . D 

Remark 8. In fact, in Proposition |7] we have that p(Z?^) = s. Again, we may assume that 
a = (0, ...,0) G N". Then, if p(Z?^) =t>s, we may suppose that Zj^* = gi, for all i= 1,.. .,f. 
We have a < b' = {b\, ...,Z7^,0, ...,0) < Z?^, p{b') = s and we get a contradiction with the 
minimality of Z?". 

Definition 9. Let a G 5<5. We define the set 

5-5(a) := {x G B>s '■ a^x, p{x) = s}. 
Theorem 10. Assume hdepthM > s. Then there exists a Hilbert partion 

^■.HM{X)^g=j^Q[a\b']{X) 
1=1 

such that if p{a') < s then b' G B-^{a). 

Proof. Since hdepthM > s, we have a partition on HM{X)^g 

^://M(X)^,= £e[a\Z7'](Z) 
1=1 

with p{b') > s. If there exists a-^ such that p{a-^) < s and b-^ is not minimal we apply 
Proposition |7] to the polynomial induced by the interval [a-^ ,bJ] and use the Remark [8] to 
end the proof. D 

Example 11. Let R = K[Xi,X2] with deg(Xi) = (1,0) and deg(X2) = (0, 1). Let M = 

R(B {Xi,X2)R. Then we may choose g = (1,1) and 

-^m(^1,^2)^(1,1) = 1+2X1+2X2 + 2X1X2. 

In order to use Corollary |5]to get that hdepthM > 1 (for details see [iT4l Example 3.5]), 
one has to compute a full Hilbert partition, for example the following 

^1 :(1 +Xi +X2 +X1X2) + (Xi +X1X2) +X2. 



AN ALGORITHM FOR COMPUTING THE MULTIGRADED HILBERT DEPTH OF A MODULE 7 

Since in this case ^ = 1, we have that 5<i = {(0,0)} and5=i((0,0)) = {(1,0),(0, 1)}. 
By using Theorem [TOl we only have to cover (0,0) with an interval ending in an element 
of 5=1 ((0,0)), and the computation is reduced to obtaining one of the following two 
possible covers: 

C:i:(l+Xi), C:2:(l+Z2). 
4. An algorithm for computing the multigraded Hilbert depth of a 

MODULE 

In this section we describe a recursive algorithm for computing the multigraded Hilbert 
depth of a module. The algorithm is presented in the form of a function which will be 
called recursively, thus realizing a backtraking search for a Hilbert partition of a given 
hdepth. The algorithm may also be used directly for computing Stanley depth in the case 
of a factor of monomial ideal. See also [[TSl Algorithm 1], for a non-recursive algorithm 
for computing Stanley depth in the case of a factor of a monomial ideal. 

Algorithm 1: Function which checks if hdepth > s recursively 

Result: true if hdepthM > s 
Boolean CheckHilbertDepth(^,5,P) ; 
begin 

ifP^N[Xi,...,X„] then return false; 
Container E =FindElementsToCover(^,5,/'); 
if size (£■) = then return true; 
else 

for /=begin(£') to?=end(£') do 

Container C[?] :=FindPossibleCovers (g,5',i',£' [/]) ; 
if size (C[i] ) = then return false; 
for j=begin(C[?']) toj=end(C[i]) do 

Polynomial P(X)=/'(X)-2[£[/],C[/][j]](X); 

if CheckHilbertDepth(g,5,/') =true then return true; 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 

11 



return /a/5e; 



At each call the function CheckHilbertDepth checks one interval of type [a, b] to see 
if the polynomial induced by it may be part of a suited Hilbert partition. All possible 
intervals are checked in a backtraking search. A node of the searching tree is represented 
by a polynomial P. Below we describe the key steps. 

• line 1 . If the polynomial P does not have natural numbers as coefficients (positive 
coefficients) then it is not a sum of polynomials induced by intervals and is not a 
node in the searching tree. 

• line 2. In this step 5<5 is computed and stored in a container. The container should 
provide some basic access functions (for example we want to ask for its size). 

• line 3. If 5<^ is empty then we are done. We have reached a good leaf of the 
searching tree. 
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• line 4,5,8. We generate and investigate all the children of the node P. 

• line 5,6. In this loop for each a E 5<i we compute the set ^^^(a) (here we use 
Theorem [Toll. 

• line 7. If ^-^(a) is empty we are in a bad node and we should go back to the 
previous node. 

• line 9,10. The child P is generated in line 9 and investigated in the recursive call 
at line 10. 

• line 11. If we have reached this point then our search in this node has failed and 
we should go back to the previous node. If we are at the root then hdepth < s. 

We conclude this section with a remark on the functions FindElementsToCover and 
FindPossibleCovers. At each node they should compute the sets 5<^ and 5=i(a) for all 
a E B^s- For a practical implementation of the algorithm it is quite inefficient to compute 
them at each node. It is likely better to adjust them for the new generated child P and pass 
them down as input data for main recursive function CheckHilbertDepth. 

5. An algorithm for computing the Stanley depth in a special case 

In this section we further assume that dim^rMa < 1 for all a E Z" and we modify Algo- 
rithm[T]for computing the Stanley depth in this case. The algorithm supplementary checks 
if the Hilbert partition computed by Algorithm [T] is inducing a Stanley decomposition. 

Algorithm 2: Function which checks if sdepth > s recursively 
Data: ^ e N", 5 G N and a polynomial P{X) = HM{X)^g E N[Xi, ...,X„] 
Result: true if sdepthM > s 
Boolean CheckStanleyDepth (g,s,P) ; 
begin 

ifP^N[Xi,...,X„] then return false; 
Container E =FindElementsToCover(^,5,i'); 
if size (£■) = then return true; 
else 

for ?=begin(£') to/=end(£') do 

Container C[i] :=FindPossibleCovers (g,s,P,E[i] ) ; 
if size (C[i\ ) = then return false; 
for j=begin(C[?']) toj=end{C[i]) do 
while a E^[E{i),C[i\[j]] do 
I if i^[Zc[,-][j]] n AnnMfl ^ then return false; 

Po\ynom\a\P{X)=P{X)-Q[E[i\,C[i\[j]]iX); 

if CheckStanleyDepth (^,5,/') =true then return true; 

return false; 

The only difference from Algorithm 1 appears at lines 2,3. Here we check if the Hilbert 
decomposition we found is a Stanley decomposition. For this we use [fT4l Proposition 
4.4]. The only thing to prove is that the conditions at lines 1, 3 assures that P is inducing 
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a Stanley decomposition. Assume that for all a E ^[E{i),C[i] [j]] we have that _^[Zc[,][^]] fl 
Ann Ma = 0. Let 7^ m^ G M^. Since dimjcMa = 1 we have that Annm^ = AnnM^, 
so ^[Z(;[,][y]] n Annm^ = 0. Then maK[Z,^^^y^] is a Stanley space. Finally, since all the 
coefficients of P are < 1, the condition at line 1 assures that they do not overlap. 
We end with a vague remark. It is easy to see that for two intervals 

[a,-, bj] n [aj, bj] 7^ <^=^ a,- V aj < bt A bj. 

Since in this particular case the intervals do not overlap, for a practical implementation 
of the algorithm one may take advantage of this fact by saving the intervals and replacing 
the test needed at line 1 . 
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